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Abstract 

We investigate the Kosterlitz-Thouless transition for hexatic order on a fluc- 
tuating spherical surface of genus zero and derive a Coulomb gas Hamiltonian 
to describe it. In the Coulomb gas Hamiltonian, charge densities arises from 
disclinations and from Gaussian curvature. There is an interaction coupling 
the difference between these two densities, whose strength is determined by 
the hexatic rigidity. We then convert it into the sine-Gordon Hamiltonian and 
find a linear coupling between a scalar field and the Gaussian curvature. Af- 
ter integrating over the shape fluctuations, we obtain the massive sine-Gordon 
Hamiltonian, which corresponds to a neutral Yukawa gas, and the interaction 
between the disclinations is screened. We find, for Ka/k 3> 1/2 where Ka 
and n are hexatic and bending rigidity, respectively, the transition is supressed 
altogether, much as the Kosterlitz-Thouless transition is supressed in an infi- 
nite 2D superconductor. If on the other hand Ka/k <C 1/2, there can be an 
effective transition. 

PACS numbers: 05.70.Jk, 68.10.-m, 87.22.Bt 



1 



Typeset using REVTf^X 



Recently, the Kosterlitz-Thouless (KT) transition for hexatic order on a free fluctuat- 
ing membrane has been investigated JIJ. A flat rigid membrane can have quasi- long-range 
(QLR) hexatic order at low temperature and undergo a KT disclination unbinding tran- 
sition to a disordered high-temperature phase. A fluctuating membrane can also have 
QLR hexatic order || at low temperature. At high temperature, a fluctuating membrane 
has no internal order and can be charaterized by a bending rigidity k. At length scales 
less than the persistence length £ p = ae 47rK//3T , where a is a molecular size and T is the 
temperature, the membrane looks flat; at longer length scales, it is crumpled. However, 
at low temperature, hexatic order stiffens the bending rigidity so that the bending rigidity 
approaches a constant times the hexatic rigidity [0]- Thus, the hexatic membrane is 
more rigid than a fluid membrane, and it is said to be crinkled rather than crumpled. A 
fluctuating hexatic membrane can undergo a KT transition from the crinkled hexatic to the 
crumpled fluid state. For fixed large k, there is a disclination melting to the crumpled fluid 
phase as temperature is increased, and at fixed Ka, there is a transition to the crumpled 
fluid phase as k is decreased. 

In this paper, we extend a study of the KT transitions to a fluctuating surface of genus 
zero. Ovrut and Thomas discussed the structure of the KT transition of a vortex-monopole 
Coulomb gas on a rigid sphere and show it is the same as in the planar case, i.e. the 
KT transition temperature on a rigid sphere is the same as that on the Euclidean plane 
; T S p£ ere = T^J ne = tt/2K a ||. We investigate the effect of thermal shape fluctuations 
of a genus zero surface on the KT transition in the limit /3k ^ 1. In this limit, we can 
parametrize the surface by its radius vector as a function of standard polar coordinates 
u = (0,0) = fi; 

R(fi) = R(l + p(ty)e r , (I) 

where e r is the radial unit vector and p(fl) measures deviation from sphericity with radius 
R. This parametrization is a "normal gauge" . To make the equations simple, we map this 
parametrization onto the unit sphere parametrization with R = 1. Later when we analyze 
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interaction between two disclinations, we will recover this length scale. Associated with 
R(O) is a metric tensor g aj g(0) = <9 Q R(0) ■ dgR(O) and a curvature tensor K a p{VL) defined 
via K al 3(Q) = N(O) ■ d a dpR(Q), where N(O) is the local unit normal to the surface. From 
the curvature tensor K a p, the mean curvature, H, and the Gaussian curvature, K, are 
defined as follows: 

H = ^K Pa , K = det g aX K X p, (2) 

where g al3 is the inverse tensor of g a p satisfying g aX g\p = 6%. 
To describe hexatic order, we construct the tangent vectors 

to = dgH, t^ = c^R, (3) 

where d a = d/du a , u = (0,<f>), and introduce orthonormal unit vectors ei and e2 at each 
point on the surface. Then, ei(d)-te(Q) = cos 0(0) defines a local bond angle 0(0). Hexatic 
order is then described by the local bond order parameter m(0) = cos 0e! + sin 0e 2 , where 
0(0) has 6-fold symmetry. Note that since 0(0) depends on the choice of orthonormal 
vectors ei and b%, any spatial derivatives for m must be covariant derivatives. 

In the continuum elastic theory, it is now well established that the long-wavelength prop- 
erties of a fluctuating hexatic membrane are described by the Helfrich-Canham Hamiltonian 
^hc @ an d the hexatic free energy Ha [§]• The Helfrich-Canham Hamiltonian can be 
expressed as a sum of three terms, 

Hhc = W K + 'Hg + ^a- (4) 
The first term is the mean curvature energy, 

= ^|rf0((V 2 + 2)p) 2 + O(p 3 ) (5) 

where y/g = det g a p, 2H = is the twice the mean curvature, Hq is the spontaneous mean 
curvature which is equal to 1 for the sphere, and the second form is valid for the normal 
gauge. The second term is the Gaussian curvature energy, 
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(6) 



where K = det Kg is the Gaussian curvature. This term is a topological invariant depending 



only on the genus of the surfaces due to the Gauss-Bonnet theorem; 



d 2 u^K = 4tt(1 - 77) = 2ttx, 



(7) 



where 77 is the number of handles and x — 2(1 — rj) is the Euler characteristic. Since we will 
consider surfaces of fixed genus, we will drop this term. Finally 



is the surface tension energy. We are mostly interested in free membranes for which the 
renormalized surface tension obtained by differentiating the total free energy T with respect 
to the total surface area A [o R = dJ-'/dA), is zero. Since there are entropic contributions 
to a R as well as contributions from internal order, the value of the bare surface tension 
a will have to be adjusted to keep a R zero. In what follows, we will ignore 7i a with the 
understanding that it is really present if we want to keep track of how a R actually becomes 
zero. 

The hexatic free energy is the contribution to the energy from fluctuations in the local 
bond order parameter. Since the hexatic order parameter m has a fixed magnitude and 
there are no external fields aligning m along a particular direction, the lowest nontrivial 
contribution to the energy associated with m arises from its gradients, 



where D a is a covariant derivative since the bond order parameter is frustrated by the 
rotation of tangent vectors that occurs under parallel transport on a curved surface. The 
amount of frustration is given by the gauge field A a , i.e. the covariant derivative of e a in 
direction a defines the gauge field A a . Under parallel transport in direction du a , each e a is 
rotated by an angle A a du a . Thus the gauge field A a is defined by 




(8) 




(9) 
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D a e a = -A a e ab e b , (10) 

where e ab is the antisymmetric tensor with e 12 = — £ 2 i = 1 an d A a e ab is called the spin- 
connection and describes how the basis vector e a rotates under parallel transport according 
to the Gaussian curvature K of the surface. In fact, A a is related to K. The curl of the 
gauge field A a is the Gaussian curvature; 

l a/3 D a Ap = K, (11) 

where the antisymmetric tensor defined via 

la/3 = N ■ (t Q x tp) = y/ge a/ 3, 

.a/3 _,aa'P/3'. flO\ 
1=9 9 la' P', (l^J 

In terms of the local bond angle G(f2), the covariant derivative of m writes 

D a m = (D a m a )e a + m a (D a e a ) 
= (D a m a )e a - m a A a e ab e b 
= (D a Q)(— sin 6e! + cos 0e 2 ) 

— A a (cos Ge 2 — sin Ge^ 
= (D a & - A a )m ± , (13) 

where m_|_ = — sin Ge! + cos Ge 2 satisfying m ■ mj_ = 0. Then the hexatic energy writes 

H A = l -K A J d 2 u^g^{d a Q - A a )(d p e - Ap). (14) 

Thus we have the Hamiltonian Ti = H K + Ha to describe fluctuating hexatic membranes. 

To gain some physical understanding of a spherical hexatic membrane, we examine the 
ground states. Since we are interested in the limit (5k >> 1 and in this limit H K dominates, 
we first minimize Ti K over the shape fluctuation field p which gives p(Q) = and then we 
minimize T~La over Q with p(f2) = and find 



6G(n) 



= -±=d b g oa \d a e° - Al) = 0, (15) 
e=e° v<? 



where the superscript stands for the rigid sphere with p(Q) = 0. In Ref. ||10|| , Lubensky 
and Prost show that in the ground state 12 disclinations of strength 2ir/Q are arranged at 
the vertices of icosahedron inscribed in the sphere. A disclination at u = u» with strength 
ft gives rise to a singular contribution G>o ing to 0° satisfying 

^du«d a er s = fi , (16) 

where T is a contour enclosing Uj. Thus, in general d a Q° = d a & + where Bq is nonsin- 
gular, = d a Q S Q ns , and 

r P D a vl = n°(Q), (17) 

where 

^(fll^pfO-fl,), (18) 

where is the disclination density in the ground state and flj's are the coordinates of the 
vertices of icosahedron. Since (d a Q° — A°) satisfies D a (d a Q° — A°) = 0, it is divergence-less 
and purely transverse. Accordingly (d a Q° — ) can be written in terms of the curl of scalar 
fields and by applying the operator ^P a D@ to (d a Q° — we find these scalar fields to be 
related to the Gaussian curvature K of the rigid sphere and the ground state disclination 
density on the rigid sphere, 

^ a D p {d a Q° - A° a ) = i Pa Dpv Q a - ^ a D p Al 

= n°-K , (19) 

where K is a Gaussian curvature of the rigid sphere and n° is the disclination density in 
the ground state. 

Now taking into account the bond angle fluctuations around 9° and the shape fluctua- 
tions around the sphere, 

e = e° + e, A a = A° a + 5A a , (20) 
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the full Hamiltonian writes Ti = Hq + 57i 



n = \k a J dQ(d a Q° - A° a )(d a Q° - A° a ) 
SH = ^kJ dQ ((V 2 + 2)p) 2 

+\k a J dQ(d a Q - 5A a )(d a Q - 5A a ) + 0(p 3 ). (21) 

The angle fluctuation field can also have disclinations of strength q = 2n(k/6) where 

k is an integer, due to the thermal fluctuation |TT|. Thus, d a Q can be decomposed into 
singular and nonsingular parts; d a <d = <9 a 6" +v a where 0" is nonsingular, v a = d a Q sms and 

-f f, D a v fi = n(n), n(n)=Y f g i s(n-ni), (22) 

i 

where n(Q) is the thermally-excited disclination density with disclinations of strength at 
flj. The vector v a can always be chosen so that it is purely transverse, i.e. D a v a = 0. In 
the hexatic Hamiltonian, d a Q always occurs in the combination (d a Q — SA a ). The spin- 
connection 5A a can and will in general have both a longitudinal and a transverse component. 
However, one can always redefine 0" to include the longitudinal part of 5A a . This amounts 
to choosing locally rotated orthonormal vectors ex(u) and e 2 (u) so that D a 5A a = 0. Thus 
we may take both v a and SA a to be transverse and the hexatic Hamiltonian 

X -K A J dtt{d a Q 11 + v a - 5A a ){d a e 11 +v a - 5A a ) 
= H\\+H ± , (23) 

can be decomposed into a regular longitudinal part, 

H\\ = \k a J dQd a Q ll d a Q 11 , (24) 

and a transverse part, 

H± = \k a J dn(v a - 5A a )(v a - 6A a ), (25) 
where the cross term / dQ(v a — 5A a )d a Q^ is dropped since D a {y a — A a ) = 0. 
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It costs an energy e c (k) to create the core of a disclination of strength k. (We assume 
for the moment that the core energies of the positive and negative disclinations are the 
same. See, however, Refs. p|,p!2| and Ref. [0.) Thus, partition sums should be weighted 



by a factor y^ = e~ l3tc ^ for each disclination of strength k. Since e c {k) ~ k 2 , we may 
at low temperature restrict our attention to configurations in which only configurations of 
strength ±1 appear. Let N± be the number of disclinations of strength ±1 and let Uj± be 
the coordinate of the core of the disclination with strength ±1 labeled by i. The hexatic 
membrane partition function can then be written as 

Z(k, K a , y) = Tr v y N j PR / pe" e -' 3m «e-^ll +w - L >, (26) 

where y — yi, and N = N + + N_. 7ij_ depends on all of the disclination coordinates fi„± 
where v = 1, 2, • ■ ■ , N±, and Tr v is the sum over all possible disclination distribution with 
the topological constraint 



where a 2 is a molecular solid angle. The Kronecker factor Sn + ,N- in Tr v imposes the topo- 
logical constraint that the total disclination strength on a sphere is 2 since with N + = AL 
we have 12 ground state disclinations with the strength 1/6 giving the total disclination 
strength 12 x (1/6) = 2. 

The hexatic model of Eq. (|26|) can easily be converted to a Coulomb gas model using 

^D a {v p - 5A P ) =n-6K = C } (28) 

which follows from Eq. (|TTD and Eq. (^) where 5K is the deviation of the Gaussian curvature 
from the rigid sphere. The quantity C = n — 5K is a "charge" density with contributions 
arising both from disclinations and from Gaussian curvature. Equation ( p8|) implies 



v a - 5A Q = --yjDp^C, (29) 

where we used 'j a \D Xr y al3 D a = —A and A = D a D a = (1/ \fg)d ay /gg af3 dp is the Laplacian 
on a surface with metric tensor g a p acting on a scalar. Recall [Eq. (@)] that ^J 3 rotates 



a vector by ir/2 so that (v a — 5A a ) is perpendicular to Dp(— A) X C and is thus manifestly 
transverse. Using Eq. in Eq. (|25|), we obtain 



where 



Z = Tr v W 2?R / P0ll e -^-^i|-^ C) (30) 

H c = ~K A J dQ^^C(Q)^-^C(Q) 

= \ka\ dQdQ'C(Q) ( ^" D ^ Dx 5{Q. - ft') J C(O') 

= ^jr A y dQdQ'C(Q) ~ fi')) (31) 

is the Coulomb Hamiltonian associated with the charge C. Since the longitudinal variable 0" 
appears only quadratically in Ti,\\, the trace over 0" can be done directly giving the Liouville 
action [15|] arising from the conformal anomaly; 



j V Q\\ e -m =e-^ L , (32) 

where 

(3H h = -— r fdQ-— I dQdQ'K{Q) (—7-6 (to - Q')) K(Qf). (33) 
8na 2 J 247T J V A / 

The Coulomb gas partition function can thus be written 

Z = Tr v y N J VRe-^-^-^c . (34) 

The Coulomb gas model can be converted following standard procedures into a sine-Gordon 
model. The first step is to carry out a Hubbard-Stratonovich transformation on j3Hc'- 

where the Liouville factor eP n ^ is needed to ensure that e~^ Wc be one when C = 0. Inserting 
this in Eq. fl3~4]), we obtain 

Z = Tr v y N J pRp$ e HW«HW* e i /<«(n-«0* ( 36 ) 



where 



= ^{PKaY 1 J dnd a $d a $. (37) 

The only dependence on disclinations is now in the term linear in n. Thus to carry out Tr v , 
we need only to evaluate 



Tt v 3/ Jv e < / dnn * 

V 1 X., „ .,*++*- ( f ^ 27r»(n)/6 ] + ( f ^ -!bri*(n)/6 

h NANA \1 a? & [J a* 6 



N+,N. 



du ( r dVt f / rdn 



V - I — I,, I — P W*(n)/6]-w} I ^ -i{27r[$(n)/6]-o;} 

^N+W-lJ 2n{ y J a* 6 ) { V J a'* ) 

_ /" du^2y/a 2 ) J dnco S [27r($/6)-w]_ (gg) 
J 27T 

Thus 

Z = J — J D$ ypRe" /3HK e _/3W *e (2y/a2) / «««»PM*/e)-u] e -i / dn<f5ic_ ^ 

We can now change variables, letting $ = (6/27r)($' + cj). The term linear in the Gaussian 
curvature then becomes 

- i J dQ5K^-(uj + $') = J dm'SK, (40) 

where we used / dQ6K = 0. The integral over uo in Eq. fl39|) is now trivial, and dropping the 
prime we obtain 

Z = j T?$ j DRe'^e^, (41) 

where 



C = l -{f3K A )~ l (^) 2 / ^<9 Q $<9 Q $ 



~ J dVtcos^-i^- / dfi<MK (42) 



is the sine-Gordon action on a fluctuating surface of genus zero. The first two terms of 
this action are the gradient and cosine energies present on a rigid sphere. The final term 
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provides the principal coupling between <3> and fluctuations in the metric. It is analogous 
to the dilaton coupling [[HJ of string theory though here the coupling constant is imaginary 
rather than real. Note that the Liouville action is not explicitly present in Eq. (|4T|). 

In the regime (3k ^> 1, we can truncate the higher order terms in p. In the normal gauge, 
the partition function becomes 

1 r .„ / o , \2 



Z = J X?pX?$exp --/3k J dQ ((V 2 + 2)p) 

--/3F / dft(V$) 2 + 4 [dQcos® 
2 J a 2 J 

J dQ$(V 2 + 2)p (43) 
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where (3Y = 36/(47t 2 (3Ka) and we used SK = (V 2 + 2) p. To lowest order in p, the shape 
fluctuation field p is linearly coupled to the scalar field $ which is the conjugate field to 
the disclinations. In Ref. M, we have shown the similar coupling in the fluctuating flat 
membrane is quadratic in the shape fluctuation field. 

Integrating over the shape fluctuation field p gives the effective Hamiltonian for the 
conjugate field to the disclinations 



Z = Jv<& exp --/ST J dQ ((V$) 2 + /i 2 $ : 



4 / dQ cos $ 



a 2 



(44) 



with p? = Ka/k. This is the massive sine-Gordon theory. The shape fluctuations induce 
the mass term for $ field and screen the Coulombic interaction between the disclinations 
giving the Yukawa interaction between them. This partition function is equivalent to that 
of the Yukawa gas Hamiltonian on the rigid sphere with radius R; 



If 1 

^Yukawa = ^P k a / dQn(Q) — ^n(Q) 
2 J — + u z ) 



where 



+ p z 

(45) 

A id 



21 + 1 

G(a-^) = E /(j + 1)+/i2 ^(cos^) 
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7T 



P v (- COS Uij), 



(46) 



cos 



where P„(-) is the Legendre polynomial with degree v = — 1/2 ± — 4/i 2 )/2 and cjjj is 
the angle between two disclinations at ilj and Qj. For < /i < 1/2, degree of Legendre 
polynomial, u, is real and the length scale introduced by A^ = (/i/i?) -1 is larger than the 
system size, 2R after recovering the original length scale by mapping the unit sphere back 
to the sphere with the radius R. On the other hand, if fi > 1/2, v is a complex number, 
v = — 1/2 ± it where r = (^Afi' 2 — l)/2 and X d < 2R. The length scale introduced by 
Xd — Rf/i may be interpreted as the Debye screening length arising from shape fluctuations. 

The interaction energy between two disclinations % and j at positions f2j and Qj with 
strength and qj is given by qiqjG(dij) where d^ = 2Rsm(iVij/2) is the chordal distance 
between two disclinations on the sphere with radius R. The interaction G(dij) has the 
following limiting forms: 



Following the analogy of the 2D Coulomb gas, when the screening length is much larger 
than the system size, A^ ^> 2R, the induced mass term arising from shape fluctuations is 
irrelavent for the KT transition and the KT transition temperature is given by T c = 
for /i 1/2. However, for /j > 1/2, the screening length is shorter than the system size, 
\d 2R, and the mass term is relavent for the KT transition and changes the universality 
class of the system. There is no KT transition at non-zero temperature. The disclinations 
are always unbound at non-zero temperature and the KT transition temperature vanishes. 
Thus we find the crossover at /i — 1/2. 

In conclusion, we present the effect of shape fluctuations on the interaction of the disclina- 
tions on a spherical surface with genus zero. We have confirmed that the screened interaction 
is of the same form as the vortex line interactions in type-II superconductors. In these su- 
perconductors, screening of vortex line interactions drives the Kosterlitz-Thouless transition 



G(dy) 



— < 



-\ In (dij/2) , fi' 1 > 2, dij < 2R 
-\ In (dij/2) , fi- 1 < 2, dij < \ d 
e~ di ^ Xd , fi- 1 < 2, d^ > X d 



(47) 
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temperature to zero for an infinite superconductor in zero magnetic field [17]]. Likewise, the 
screening of the disclination interaction on the fluctuating spherical surface drives the KT 
transition temperature to zero for /i -1 <C 2 in which the screening length is shorter than the 
system size. However, when ^> 2, the effect of shape fluctuations is irrelavent and the 
effective KT transition occurs at the finite temperature. 
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